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1 Introduction 



The dynamics of a system of strongly correlated electrons on a lattice is most directly 
described by a hamiltonian that acts on the microscopic quantum states of the system. 

The simplest example is the hamiltonian of the Heisenberg XXX magnet, which 
describes spin-spin interactions of electrons that are nearest neighbors on the lattice. 
In this model the electrons are not mobile; one assumes that each lattice site has 
either a spin-up or a spin-down electron, so that the total number of independent 
states, for a lattice of L sites, is equal to 2^. 

If one is interested in modelling phenomena such as (high-Tc) superconductivity, 
one allows empty or doubly occupied sites and includes in the hamiltonian kinetic (or 
hopping) terms which allow the electrons to move. A simple example is the so-called 
t-J hamiltonian, which acts in a space of states (of dimension 3^) in which doubly 
occupied sites are forbidden. Besides the kinetic terms, the t-J hamiltonian contains 
nearest neighbor density- density and spin-spin couplings. 

A further refinement can be made by allowing both empty and doubly occupied 
sites (enlarging the basis of the Hilbert space to 4^ states). The probability of finding 
doubly occcupied sites can be adjusted through an on-site potential term, the so- 
called Hubbard potential, which contributes an amount U of potential energy for each 
doubly occupied site. In the so-called Hubbard model there are no interactions besides 
this potential term, but it is very natural* to include nearest neighbor interactions 
such as the density-density (V^) and spin-spin (J) couplings of the t-J model and 
additional interactions such as a bond-charge term {X) and a pair-hopping term [Y). 
One certainly expects that a hamiltonian that has each of these couplings tuned to a 
realistic strength will give an accurate description of the physics of the electron gas. 

In formula, this general hamiltonian takes the following form 

<ij>,a 

+V Y: {n^ - l){n, - 1) + J E iS^S] + ISjS, + ^SjS,) 

<ij> <ij> 
<ij>,cr 

*See for example Hubbard's original paper 0|, where estimates for the strengths of some additional 
interaction terms are given. 
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+^ J2 iwl-lCj,-lCj,l + cJ^iC}_„iCi,_iCi,i) 



+f/EKi-|)K-i-i) ' (1-1) 

i 

where < ij > denotes the nearest neighbor relation on the lattice. The operators 
cl^/ci^a- {i = 1,2, cr = ±1) create/annihilate electron states. They satisfy 

{c|,y, Cj^r] = Si,jSa,T- The space of states is a Fock space with the pseudo vacuum 
I 0) defined by Cj^o- | 0) = 0. The operator rij o- = c|^Cj^o- is the number operator for 
electrons with spin a on site i and rii = n^^i + nj _i. The SU{2) spin operators are 
given by = - ra^ _i), Sj = 4_iCi,i and Si = cJiQ _i. 

The general hamiltonian ( |1 . 1| ) is obviously hard to analyze. However, life becomes 
simpler if one assumes certain relations between the coefficients of the interaction 
terms. In general, exact results are most easily obtained in models that have a high 
degree of symmetry, and one may therefore look for special cases where the hamilto- 
nian (|1 . 1|) admits symmetries that go beyond the SU{2) spin rotations. Following this 
line of thought, we observe that for each choice of the kinematical constraint on the 
electron states (allowing 2,3 or 4 states per site), there exists a hamiltonian that is 
maximally symmetric, i.e. that is invariant under all global, unitary rotations among 
the elementary electron states. The Lie (super-) algebras that describe these symme- 
tries are SU{2), SU{2\1) and SU{2\2), respectively, where the grading expresses the 
fact that the single electron states have fermionic statistics. 

The first two maximally symmetric models are well-known: the SU{2) symmetric 
case is the XXX Heisenberg magnet and the model with SU (2|1) symmetry is the t-J 
model at the supersymmetric point J = 2t. The model with SU {2\2) symmetry, which 
is the logical third in this list, has been considered only recently 0]. Its hamiltonian 
is of the form ( |1.1| ) where V = —t/2, J = 2t, X = t, Y = —t, and U = —dt, where 
d is the number of nearest neighbors of each lattice site. The physical properties of 
this model are quite remarkable: it was shown in ^ that the groundstate in 1,2 or 3 
dimensions is superconducting, even when a moderately repulsive Hubbard potential 
term is added to the hamiltonian. 

The superconductivity in the SU{2\2) model comes about through the Bose- 
Einstein (BE) condensation of local electron pairs (bipolarons), which is in contrast to 
the usual BCS mechanism. It has been observed (see for example 0, §) that the prop- 
erties of 'BE superconductors' agree rather nicely with the observed phenomenology 
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of the high-Tc superconductors (this is for both the superconducting and the normal 
state). We therefore expect the SU{2\2) model, which can be viewed as an 'ideal BE 
superconductor'^, can play an important role in the quest for a theory for high-Tc 
superconductivity. We refer to a forthcoming paper for further discussion of the 
physics of the model. 

From a mathematical point of view the three maximally symmetric electron models 
are relatively easy to work with: in particular, they share the property of being exactly 
solvable in 1 dimension. (The Hubbard model, which is not maximally symmetric, 
also has this property The exact solution of the SU {2\2) model employs a three- 
fold nested, graded Bethe Ansatz. It was presented in [|], see also and has been 
further worked out in |1T0|. 

The Bethe Ansatz being an Ansatz, there is of course the issue whether or not it 
leads to a complete solution of a given model. A more precise question is whether the 
Bethe Ansatz, together with the global symmetries of the model, provides a complete 
set of eigenstates of the hamiltonian (complete in the sense that these states form 
a basis for the Hilbert space of the model). For the XXX model this issue was 
already discussed in Bethe's original paper For the Hubbard model and for the 
supersymmetric t-J model completeness has been established only recently |T2|, |13[. In 
this paper we shall extend these results to the SU{2\2) extended model and explicitly 
discuss the counting of the eigenstates. 

The paper is organized as follows. In the next section we briefly review the Bethe 
Ansatz and we discuss the so-called string solutions to the Bethe equations. We derive 
a formula for the number of solutions that can be constructed by combining a specific 
set of string solutions, and we specify the dimensions of the SU{2\2) multiplets that 
are present in the spectrum. The counting of the eigenstates is discussed in detail in 
section 3. 



2 Analyzing the Bethe equations 

The exact solution of the SU {2\2) extended Hubbard model on a one- dimensional 
lattice starts from the observation that the hamiltonian Hq (which is as in ( |1 . 1| ) 

^ 'Ideal' in the sense that momentum-zero electron pairs only interact via a chemical potential 
term. For non-zero momentum pairs this is no longer the case. 
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with t = 1, V = —1/2, J = 2, X = 1, Y = —1 and U = —2) acts as minus the sum 
of graded permutations 11? • of the electronic states on neighboring sites i, j : 

^0 = - E % ■ (2-1) 

<tj> 

Lattice models for which the hamiltonian acts as a (graded) permutation were first 
considered in |T3]. We shall here use the so-called FFBB representation of the Bethe 



Ansatz solution, which was presented in detail in |^ (see also 0). The details of the 
Bethe Ansatz approach are rather complicated, but the final result is easily explained 
as follows. The (nested) Bethe Ansatz provides a multi-parameter family of states, 
written as 

which turn into exact eigenstates of the hamiltonian Ho when the spectral parameters 



{Afc, Xf\ X^} solve the following set of algebraic equations 



M + h) j=i Afc - x'-p - ^ i=i,i^k 



j = l Ak — Aj — 2 l=l,l^k 

Ni+Ne \ _ , i \(2) _ , i 

-pj- ^ 2 _ TT "^m -^3^2 ■ _ i 

11 . .(1) , 11 ,(2) ,(1) i ' J 

k=l ^k ~ '^j ~ 2 ™=1 ^"^ i 2 

n =n ^s, m-i.-.-iv.. (2.3) 

i=l,i^m — Am — I j=l Am ~ +2 

These equations are called the Bethe equations for this specific Bethe Ansatz. By 
the symbols A^;, A^e, and we denote the number of local electron pairs, of single 
electrons and of single spin-down electrons, respectively. The total numbers of spin-up 
(N-^) and spin-down (A^|) electrons are thus given by 

N^. = Ni + Ne-Nd, N^ = Ni + Nd. (2.4) 

The energy of the eigenstate ( |2.2| ) for a given solution {Xk, Xf\ X^^^} of the Bethe 
equations is given by 

E= Y: tt-T-L. (2.5) 

k=l ^fc + 4 
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Before we come to the description of the solutions to the Bethe equations, we 
wish to make two important remarks, which are both relevant for the counting of 
eigenstates that we are going to perform. 

The first remark concerns an especially nice property of the Bethe Ansatz eigen- 
states, which is that they are lowest weight states of the full symmetry algebra SU (2|2) 
1^. This implies that each Bethe Ansatz eigenstate comes with an entire multiplet of 
eigenstates of the hamiltonian. The dimensions of these multiplets, which are derived 
by using the representation theory of the algebra SU {2\2), are given below in the 
equation ( [2.1(J| ). 

The second remark concerns the norm of the Bethe Ansatz eigenstate for a given 
solution {Afc, A^-^"*, A^-*}. For the simpler case of the non-nested bosonic Bethe Ansatz 
there exists an elegant formula for the norm of all Bethe Ansatz eigenstates JTSf , which 
in particular implies that this norm is non- vanishing for so-called regular solutions to 
the Bethe equations. However, for the present situation a closed expression for the 
norm of the Bethe Ansatz eigenstates is not known, and we have to rely on other 
arguments to find out which solutions to the Bethe equations ( |2.3| ) give rise to non- 
vanishing eigenstates. We propose the term spurious solution, for a solution to the 
Bethe equations that corresponds to a vanishing eigenstate. 

Let us now focus on the solutions of the equations (p73|). We will interpret general 
solutions of ( p.3|) as combinations of so-called string configurations, which are config- 
urations in the complex A-plane that are close to the string solutions that solve the 
Bethe equations in the limit where L tends to infinity. For the present model, the 



relevant string solutions are |TD 



A-strings : A^^' = A™ + ^(m + 1 - 2j) , j = 1, . . . ,m 
real A^^) : A^'^ 

A(i)_A(2) strings : Af , A^^' = Af + ^(3 - 2j) , j = 1,2 

A(2)-strings : A^,')'"^' = A^,')"^ + ^(m + 1 - 2j) , j = l,...,m, (2.6) 

where the string centers A™, A^"^'', A^^-* and X^^^"^ are all real numbers. Notice that the 
XW-Xi"^) strings each consist of one real rapidity A^^^^ and two complex rapidities Xi^^^ , 
J = 1,2. 
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We should now determine in what ways these string configurations can be com- 
bined into true solutions of the equations (p.3|) . This can be done by substituting the 
string solutions (p^) into the Bethe equations and extracting a set of equations for 
the string-centers. We shall focus on a solution that contains Mm A-strings of length 
m, Ni real A*-^-*, N12 A'-^-'-A*-^-' strings and K^, A -strings of length m. 

Imposing periodic boundary conditions implies that the string centers have to be 
chosen from a discrete set. We will make the (standard) 'assumption' (see below for 
the quotation marks) that this set is parametrized by a set of non-repeating strictly 
positive integers, one for each string configuration that is present in the solution. We 
will denote the integers corresponding to the string centers A™, A^^"*, and X'^^^^ by 
I^, J/3, K^, and NJJ^, respectively. The allowed ranges of these integers were studied 
in |jTO|. Following this work (with one modification, see below), we claim that the non- 
spurious solutions to the Bethe equations are parametrized by non-repeating integers 
satisfying 



0<I^<L- Y^itmnMn) + iVi + iVi2(2 - 5, 



ml ) 



n=l 



0<J/3<J2(Mrn- Km) - 



12 



m=l 



< < M1 + 2J2 Mm- Ni- N: 



12 



m=2 



< iV™ < iVi - ^ (tmnKn) , (2.7) 
m=l 

where tmn = 2 Min(m, n) — 6mn- This result implies that the total number of solutions 
with string multiplicities {Mm, Ni, N12, Km} is given by 

m=l \ -^m 



n{{Mm},N,,N,2,{Km})= n 



X 



X 



En=liMn - K^) - Nu - I 



Mi+2E^=2^n-^l-^12 



7 



X n I ^^=^^'""^"^ I . (2.8) 



m=l 



Before we continue, several comments on the validity of this formula are in order. 

A first comment concerns the 'assumption' that the solutions to the Bethe equa- 
tions ( |2.3| ) are in 1-1 correspondence with sets of non-repeating integers that satisfy 
the constraints ( p.?] ). This assumption is similar to an assumption made in Bethe's 
original work on the Heisenberg XXX chain [|1T|, which was followed in later work by 



many. However, Bethe immediately pointed out that this assumption gives a slightly 
idealized version of the situation, since the actual string content of solutions can be dif- 
ferent from the one implied by this counting. This phenomenon was carefully studied 
in ref. [|^, where it was argued that the assumption is correct up to a redistribution 



of solutions of one type into another. (For example: it was found that for increasing 
lattice length certain complex 2-magnon string solutions in the XXX model turn 
into pairs of real solutions.) With this in mind, we will continue our discussion on the 
basis of the equation ( p.7|) . 

A second comment concerns the upper limits on the integers J™, Jp, and X^. 
Our claim is that the ranges specified in ( p.7|) are such that they precisely select so- 
lutions to the Bethe equations that correspond to non-vanishing Bethe Ansatz eigen- 
states. We have tested this claim by explicitly working out the SU {2\2) quantum 
numbers of all lowest weight states for the model on L = 2, 3, . . . , 9 sites. The mul- 
tiplicities obtained are precisely reproduced by the expression (p!8|). There can thus 
be little doubt that the proposed expression is correct (the check at L = 9 involves a 
total of 262,144 states, 2,578 of which are lowest weight states). 

If one is interested in obtaining solutions to the Bethe equations (without insist- 
ing that the corresponding eigenstates be non- vanishing) , one finds a slightly larger 
domain for the allowed integers /™, J/3, and N^. This observation explains the 
difference between the inequalities ( P?7| ) and those reported in []TU[. The following set 



of rapidities provides an example of a spurious solution to the Bethe equations (|2T 
for L = 3 

Ai = ^v^, A2 = -^v^, A« = ^v^, X? = -'-V3, A(2) = 0. (2.9) 
o D 3 3 

It can be viewed as a combination of two real A's and one A''^^-A''^^ string, with Mi = 2, 
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A''i2 = 1. It is a perfectly sound solution to the Bethe equations and as such its 



existence was predicted by the analysis of [|^. However, it can easily be seen that 
this solution does not correspond to a non-vanishing eigenstate of the hamiltonian. 
The reasoning is as follows. We know that the state has quantum numbers A*"; = 0, 
= 1, = 1, and that its energy is E = 3. It is easily checked that the spectrum 
of the L = 3 model has only one state with these characteristics, and that this state 
is not a lowest weight state for SU {2\2). Since the theorem of guarantees that all 
Bethe Ansatz states are lowest weight states, we have to conclude that the state with 
spectral parameters (|2.!j|) is actually vanishing. This result has been confirmed by 



explicit evaluation of the Bethe Ansatz eigenstate JT^. Recently, the solution (|2] 



has been generalized to a large class of spurious solutions fl^ . 

Before closing this section, we come back to the multiplet structure of the space 
of eigenstates of the hamiltonian (|2.1|) . We already mentioned that the lowest weight 
theorem for Bethe Ansatz eigenstates implies that every Bethe Ansatz eigenstate 
gives rise to a multiplet of eigenstates which form a representation of SU{2\2). The 
dimensions of these multiplets can be obtained from the representation theory of 



the superalgebra SU{2\2), as explained for example in [0. It turns out that the 
dimensions can be expressed in terms of the numbers and N^. Depending on the 
value of Ni, the representations can be typical or atypical. The precise result is 

typical: dim^v^^TV; = 16(L - A^^ - A^^ + 1){N^ - A^^ + 1) 

if N^>Ni>2, N^ + Ni<L , 
atypical: dim7v^,7v^ = (8A^| + 4)L - 8A^^^ - 8A^| 



dimAT^^AT^ = 4L 



if Ni = l, 1 < A^i < L - 1 



if A^^ = A^^ = . (2.10) 



The numbers and Ni are expressed in the multiplicities {Mm, Ni, N12, Km} ac- 
cording to 

AT^ = J2 ^(Mrn - Km) - N12 
m=l 
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= £m(M^ + i^J-iVi-iVi2 . (2.11) 



m=l 



3 Counting eigenstates 

With the results of section 2, we are now ready to discuss the completeness of the 
Bethe Ansatz solution. This issue was already addressed in [|10|, where it was shown 



that in the thermodynamic limit the leading large-temperature behavior of the free 
energy is given hj F/L ~ — Tln(4). This relation indicates that the number of 
eigenstates generated by the Bethe Ansatz is of the order of 4^, but due to possible 
sub leading corrections this is not a rigorous result. 

Using ( |2.8| ) and ( [^.1U| ), we obtain the following exact expression for the total 
number of non-vanishing eigenstates that can be obtained by combining the Bethe 
Ansatz with the 5'L'^(2|2) multiplet structure 

# (eigenstates) = n{Mm, Ni, N^, Km) dimN^^N^ ■ (3.1) 

Mm,Nl,Nl2,Km>0 

In this section we will prove that for general L this sum equals 4^. If we then use 
(as is easily shown) that the eigenstates in (|3.1| ) are all linearly independent, we may 
conclude that the Bethe Ansatz does give a complete solution of the model. 

As an aside, we remark that the SU {2\2) multiplet structure of the space of eigen- 
states of Hq is a consequence of the fact that we imposed periodic boundary conditions 
on the model. If we would instead impose twisted boundary conditions, which in gen- 
eral break the SU {2\2) symmetry, we would lose the multiplet structure and obtain a 
larger set of Bethe Ansatz eigenstates (compare with [|I^). 

As a first step in the evaluation of ( |3.1D , we shall identify the string multiplici- 
ties {Mm, Ni, Ni2, Km} that lead to atypical representations. It is easily shown that 
these only arise if all multiplicities are zero (this gives the vacuum representation of 
dimension 4L) and for the choice Mi = n, Ni = n — 1, with the other multiplicities 
vanishing. The contribution of these atypical multiplets to the summation in ( |3.1|) is 
easily evaluated 



#(atypical) = 4L + V (8n + 4)L - Sn^ - 8n 
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= 2^(L2-L + 2) 



(3.2) 



Notice that for L = 1, 2, 3 this number is actually equal to 4^; for these small lattices 
all eigenstates of the hamiltonian belong to atypical multiplets. 

This was the easy part. The hard part is to determine the contribution from the 
typical multiplets. It is given by 



#(typical) = E n 



X 



X 



Mn,Ni,N-,2,Kmm=l \ Mm 



E^=i(M„-X„)-7V,2-l 



M, + 2En=2Mn-m-Nu-2 



m=l 



X 16 



L + 1-2 Y^{mMm) + ATi + 2N^2 



(3.3) 



where the * on the summation symbol indicates that we should exclude the terms 
that correspond to atypical multiplets. Our task is now to evaluate this summation 
in closed form. 

We will tackle this problem in a number of steps, where we will be using results 
that were obtained in the analogous counting problems for simpler models, notably 
the spin- 1 XXX quantum chain. 

Step 1. 

We recall the following identity, which was first written by Bethe in his original paper 
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on the solution of the spin-i XXX quantum chain, and which is vahd for < [L/2] 



Ni,N2,...=0 [[ 



M L- N+1 \ n \ n-1 



(3.4) 

Using this identity we can directly perform the summation over Km, ttl = 1,2, . . ., 
with the restrictions J2m=ii''^^rn) = ^2 and J2m=i^rn = K. If we shift the integer 
label of the Mm by one unit and use the fact that tm+i,n+i = 2 + tm,n, we can use the 
same identity (|3.4|) to perform the summation over M2, M3, . . . with the restrictions 
J2m=2{''^^m) = M and J2m=2 = 1- The resulting expression is 



# (typical) = 

Mi,M,q,Ni,Ni2,K,N2 



L-Mi-2q + N^ + N^2 
Ml 



X 



X 



X 



L - 2Mi - M -q + Ni + 2X^2 + l\ / M -q-1 
q ) \ q-1 

Mi + q-Xi2-K -l\ / Mi + 2q- Xi- Xu-2 
iVi j I Nu 

Xi-X2 + l\ I X2-I 
K I \ K-l 



{L - 2Mi -2M + Xi + 2Xi2 + 1)^ (iVi - 2X2 + 1)^ 
^ L- 2Mi -M -q + Xi + 2Xi2 + 1 A^i - iVs + 1 ' ^ ' ^ 



Step 2. 



The summation over A^2 can be performed as follows 
W^] ( Xi-X2 + l\ ( X2-l \ (iVi - 2X2 + 1) 

N2=0 



K \ K-l X1-X2 + I 



12 



(N, - 2N2 + 1) 

N2=0 
N2=0 



Ni + 1 
2K + 1 

where we used the relations 



iVi - 7V2 + 1 \ N2-I 
K j\ K-1 

Ni-N2\ I N2 
K-1 j \ K 



i^ol K-1 [ K 



^1+1 / N1-N2 

Y{N2-K)' 

N2=0 



Ni + 1 
2K 



N2 
K 



(K + l) 



K-1 

which are easily proved. The sum over K is performed as 



K-1 



Ni + 1 
2K+1 



N2 
K 



(3.6) 



(3.7) 



E 

K 



Ni + 1 \ X-K 



2K + 1 



2X-N1 + I 



(3.8) 



with X = Ml + q — N12 — 1. We then change variables from {Mi, M, Ni, q, N12} to 
{A,B,C,q,Ni2}, where 

A^2Mi + 2q-Ni- 2Ni2 , B^M-q, C ^ Mi + q- Ni- N^-l ■ (3.9) 



With that, the summation over A^i2 becomes elementary and we obtain 



# (typical) = J2 

A,B,C,q 



X 



L-2 

A-2 



L-A-B+1 



A-1 
2C+1 



16 



B - 1 
q-1 

{L-A-2B + 1] 
L-A-B+1 



(3.10) 
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Step 3. 



Before we continue, let us specify the range of the variables q and C. The summand 
in the above is nonvanishing when q > and C > 0. However, we should exclude 
the terms that give rise to atypical representations of SU{2\2); it can easily be seen 
that those correspond precisely to the choice q = C = 0. The summations should 
thus be carried out over the regions q > 1, C > and g = 0, C > 1. Performing the 
summation over C we obtain 



#(typical) = ^ 

A,B,q 



L-A-B+1 



L-2 
A-2 



xl6 



{L-A-2B+1) 



L-A-B+1 
The summation over q then gives 

#(typical) = I 

where 



B - 1 
q-1 

2^-2 if g > 1 

2^-2 - (A - 1) if g = . 

-II 



I = ^16 



{L-A-2B + lf / L-A 
L-A-B+1 [ B 



L-2 
A-2 



^16{L- A-2B + 1] 

A,B 



L-A 
B 



L-A 
B-1 



L-2 
A-2 



and 



77 = ^16(L-A + 1)(A 

A 



-11 



L-2 
A-2 



(3.11) 



(3.12) 



iA-2 



(3.13) 



(3.14) 



We can now use the following lemma, which originates from the counting of states in 
the spin-i XXX-model 



[N/2] r 

E 

M=0 - 



N 
M 



N 
M-1 



(iV - 2 M + 1) = 2^ , 



(3.15) 



14 



to obtain 



/ = 5^16 1^ ^12^-^2^-2^22^ = 4^. (3.16) 



A 



A-2 



The second term can be evaluated as 

// = ^16(L-A + 1)(A-1) 

A 

= 2^(L2-L + 2) (3.17) 



L-2 
A-2 



and we find 



#(typical) = 4^ - 2^(L2 - L + 2) . (3.18) 



Comparing with (p.2| ), we see that the contribution from the atypical multiplets pre- 
cisely cancels the second term in ( p.l8| ) so that the full result of the summation ( |3.1| ) 



is given by 

#(eigenstates) = 4^ . (3.19) 

This is the desired result. 
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